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Abstract 

The present manuscript summarizes the modern view on the problem of the graphene-metal 
interaction. Presently, the close-packed surfaces of d metals are used as templates for the prepara- 
tion of highly-ordered graphene layers. Different classifications can be introduced for these systems: 
graphene on lattice-matched and graphene on lattice-mismatched surfaces where the interaction 
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with the metallic substrate can be either "strong" or "weak" . Here these classifications, with the 
focus on the specific features in the electronic structure in all cases, are considered on the basis of 

£ 

large amounts of experimental and theoretical data, summarized and discussed. The perspectives 
of the graphene-metal interface in fundamental and applied physics and chemistry are pointed out. 
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I. INTRODUCTION 



The physics and chemistry of graphene overlayers on metallic surfaces have become the 
focus of worldwide research since the extraordinary properties of the single layer graphene 
(Fig. [l| were demonstrated in 2004 p~H3]. This topic is interesting from several points of 
view. From the practical side, as was recently demonstrated, the metallic substrate can 
be used for the preparation of graphene layers of different thicknesses with the extraordi- 
nary quality that can be transferred onto an insulating or polymer support [IH7]. In the 
latter case the obtained material was successfully used for the fabrication of flexible touch 
screens [7j. These facts together with the relative easiness of the preparation procedure 
of graphene on metals and its transfer onto the polymer or insulating support make this 
method a most promising one for the graphene industry. Additionally, in any of the techni- 
cal applications the metallic contacts to a graphene layer determine the doping of graphene 
(and, as a consequence, its transport properties) that points the importance of studies of 
the graphene /metal interface jSHTU]. 

From the other, more fundamental, point of view, the nature of bonding itself at the 
graphene-metal interface is still far from a complete understanding. For example, the weak- 
ening of the bonding between graphene and d-metals from the 4th and 5th periods of the 
periodic table can be related to the filling of the corresponding d shell. This leads to a 
decrease of the degree of hybridization between d states of the metal and the ir states of 
graphene (effect of hybridization has to fulfill two conditions simultaneously: spatial and 
energetic overlapping of orbitals). The representative example here is the difference of the 
graphene behavior on the Ni(lll) and Cu(lll) surfaces. Both interfaces are lattice-matched. 
The graphene/Ni(lll) interface is the representative example of the so-called "strong" in- 
teraction between graphene and the metallic substrate with the significant hybridization of 
the graphene ir and Ni 3d states that leads to the complete destruction of the graphene Dirac 
cone [9"HT2"]. In the case of "weakly" bonded graphene/Cu(lll) the theoretical description of 
this system [the small lattice-mismatch of 3.9% between graphene and Cu(lll) is omitted] 
shows that the graphene-derived Dirac cone is preserved in the electronic structure [9] and 
its position is defined by the doping level. These observations are explained by the absence 
of hybridization around the Fermi level (Ep) between Cu 3d and graphene it states. A similar 
approach can be used for the description of bonding mechanisms of graphene on surfaces of 
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FIG. 1: (a) Crystal structure of the graphene layer where carbon atoms are arranged in a 
honeycomb lattice. The unit cell of graphene with lattice constant a has two carbon atoms per 
unit cell, A and B. (b) Electronic dispersion of ir and tt* states in the honeycomb lattice of free- 
standing graphene obtained in the framework of tight-binding approach. These branches have linear 
dispersion in the vicinity of the K points of the Brillouin zone of graphene. (c) Band structure of 
free-standing graphene as obtained by means of DFT (a, tt, and tt* bands are marked). 

4d and 5d transition metals. 

The situation becomes more intriguing if one now considers the interaction of graphene 
with the d-metals from the same group. Good examples are the Co-Rh-Ir or Ni-Pd-Pt se- 
quences. In both of them the strength of interaction of graphene with a metallic surface is 
decreased. Usually, as a measure of the strength of the graphene interaction with metallic 
surfaces one uses the bonding energy per carbon atom, the state of the Dirac cone or the 
doping level which is defined by the position of the Fermi level with respect to the Dirac 
point. It is interesting to note that in all graphene-metal systems the bonding energy of 
graphene is too small compared to the value which can be estimated from the temperature 




FIG. 2: (a) Crystallographic structure of the graphene/Rh(lll) interface. The high-symmetry 
stacking positions of carbon atoms on Rh(lll) are marked by different symbols (see text for details), 
(b) Height variation of carbon atoms in the graphene/Rh(lll) system. The zero level is taken for 
the BRIDGE-position. 

needed for the graphene synthesis, e. g. by chemical- vapor deposition (CVD) method. Here 
the synthesis temperature is varied from 450°C for Co(0001) to 1200°C for Ir(lll), that is 
opposite to the current description [T314T5] . Following the simple estimation for the energy 
E = ksT, graphene has to be bonded stronger to Ir(lll) compared to Co(0001), that is op- 
posite to results obtained from the theoretical calculations: the graphene layer on Co(0001) 
is considered as a "strongly" interacting system with a bonding energy of 0.132eV/C-atom, 
a high doping level of graphene, and a fully destroyed Dirac cone due to the strong hybridiza- 
tion of Co3d and graphene 7r states [16j . Graphene on Ir(lll) is a "weakly" bonded system 
with a bonding energy of 0.05 eV/C-atom, very low doping level of graphene, and an intact 
Dirac cone [UJ . From the other side, taking into account the higher catalytic activity of the 
Co surface compared to Ir one we can expect that graphene will grow faster on Co(0001) at 
higher temperatures. However, this is not the case as was demonstrated in Ref. where 
the decreasing of the graphene growth rate was observed at temperatures higher than 800° C. 
Most probably this is due to the instability of graphene on Ni(lll) and Co(0001) at such 
high temperatures. At the same time the mismatch between a graphene layer and the metal- 
lic substrate can not be considered as an argument explaining the bonding of graphene on 
metals: Rh(lll) and Ir(lll) substrates both do not match to the graphene plane. However, 
graphene/Rh(lll) is a "strongly" bonded system opposite to graphene/Ir(lll) [18, 19J. 
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Another interesting question related to the graphene-metal interface is the existense or 
absence of the energy gap in the energy dispersion for the graphene ir states around the Dirac 
point. Generally, the appearance of such an energy gap is connected with the breaking of 
the sublattice symmetry in the graphene layer. However the carbon atoms of the graphene 
unit cell occupy nonequivalent positions in all cases of the graphene adsorption on a metallic 
substrate. Thus, one should expect violation of the sublattice symmetry and, as a conse- 
quence, the appearance of an energy gap in the electronic spectrum of graphene. In fact, the 
situation is not so unequivocal. Whereas, for the "strongly" bonded graphene on metals the 
appearance of such a gap is more or less obvious (due to the strong hybridization between 
valence band electronic states of graphene and metal and high level of doping) and was ob- 
served experimentally (and confirmed theoretically) for graphene on Ni(lll), Co(0001), or 
Ru(0001) [201422] . for the "weakly" bonded graphene the experimental data are inconsistent. 
For example, some experimental results for graphene on Au(lll) or Cu(lll) indicates the 
existence of such a gap [231121], whereas for the graphene/Al(lll) system experimental data 
demonstrate no gap [HI I2S] • The theoretical calculations predict the gap-less situation in all 
cases of the "weakly" bonded graphene on metals [HI [26] . 

The present article is intended to give a timely account of recent developments in the 
investigation of the graphene-metal systems, in particular the aspects of the interaction at 
the interface will be highlighted. These points will be reviewed taking into account the 
massive set of theoretical and experimental data. The paper is organized as follows. The 
introductory part contains the main questions which arise in the course of the study of 
the graphene-metal interface. The modern preparation methods as well as computational 
and experimental approaches are briefly reviewed in the second section. Sections III and IV 
demonstrate several representative examples of the graphene-metal systems which are used 
for the critical discussion of the existing problems rising in the introductory part and others. 
Conclusions are given in Section V. 

II. EXPERIMENTAL AND THEORETICAL APPROACHES 

Many references on the particular experimental and theoretical works as well as on the 
preparation methods of graphene on metals can be found in recent review articles [T3l ~ 
[15]. Here we briefly overview the main experimental and theoretical methods used for the 
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FIG. 3: STM images of a graphene layer on Ni(lll), Rh(lll), and Ir(lll). The corresponding 
LEED images are shown as insets. The capital letters in the middle STM image denote the high- 
symmetry adsorption sites of carbon atoms on Rh(lll). 

investigation of the graphene /metal systems. 

Presently there are two commonly used ways of the preparation of graphene on metals [D- 
[7j. Both of them include the reaction of the metallic surfaces with hydrocarbons. In the 
first case the metallic substrate is kept at high temperature and pressure of hydrocarbons 
during long time that leads to the dissolving of carbon atoms in the bulk metal. The cooling 
of the metallic substrate leads to the segregation of the carbon atoms at the surface of the 
metal. The careful control of temperature, pressure and the rate of cooling can produce 
high quality graphene, from one layer to multilayers. This method can be successfully used 
for the preparation of graphene layers of controllable thickness on polycrystalline metallic 
support (mainly on Ni, but was also demonstrated for Cu) and it was demonstrated that 
graphene obtained in such a way and transferred on the polymer insulating support can 
be used in the industrial applications for fabrication of gas sensors or touch screens. The 
second method uses the high catalytic activity of low-index metallic substrates and that the 
single graphene-layer growth is a self-limiting process due to the inertness of the graphene 
surface [27J. Here the surface of metal heated to the particular temperature is exposed to 
the hydrocarbon atoms at the pressure of — 10~ 5 mbar allowing the formation of single 
graphene layers of very high quality. 

The graphene /metal system is an ideal object for the surface science studies due to the 
two-dimensional geometry of the system and consequently the character of the electronic 
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states. Initially the problem of the graphene layer (or monolayer of graphite in former time) 
on metal was connected with the catalytic properties of metallic surfaces and the passivation 
properties of graphene. It was found that adsorption or segregation of carbon layers leads 
to the quenching of the catalytic activity of metallic surfaces. Correspondingly, focusing 
on the applied research, most of the studies were performed with the application of low- 
energy electron diffraction (LEED) and Auger electron spectroscopy (AES) methods with 
rare applications of photoelectron spectroscopy. Several applications of (high-resolution) 
electron-energy loss spectroscopy methods to graphene/metal systems were published with 
the aim of determination of the strength of interaction between carbon layer and metallic 
surface. With the discovery of the fascinating properties of the free-standing graphene, the 
main focus of the surface science studies was moved to the investigation of crystallographic 
structure and electronic properties of the graphene/metal systems by means of scanning 
tunneling microscopy (STM) and angle-resolved photoelectron spectroscopy (ARPES), re- 
spectively. Here, due to the fact that in many cases of the prepared graphene/metal interface 
the dispersion of graphene-derived 7r states is strongly disturbed, compared to free-standing 
graphene, a large number of works were devoted to the modification of such systems via ad- 
sorption or intercalation of different materials on top or underneath of the graphene layer. 
Here we would like to address the readers to the several recent reviews on graphene/metals 
systems where many experimental works are cited [T3] - [T5] . The corresponding references on 
particular works will be introduced later in the discussion. 

The most widely used theoretical approach, when studying the graphene/metal systems, 
is density functional theory (DFT), ranging from the X a method [28] to the recently revisited 
random phase approximation (RPA) [291 130] . Initially, the lattice-matched graphene/metal 
systems were in the focus of these studies [TT] . Parallel execution of DFT-codes, like VASP, 
available for many users allows to study the graphene/metal lattice mismatched systems, 
where graphene and metal layer form moire structures with periodicity of ca. 3 nm and the 
model unit cells contain few hundreds of atoms [31]. Still, most of the studies were performed 
employing the local density approximation (LDA), yielding reasonable binding energies (due 
to error cancellation) contrary to the standard generalized gradient approximation (GGA). 
In the latter case, the weakly bonded structures are usually predicted to be less stable 
than in reality, and sometimes they cannot even be localized as a stationary point on the 
potential energy surface [2H1 EO] • (It is necessary to note here, that density functional theory 
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itself is capable of providing the exact solution to the Schrodinger equation, including long- 
range correlations - the dispersion. The root of the rather unsatisfactory description of the 
intermolecular interactions in DFT lies in the approximations made in the DFT functionals 
of all types: LDA, GGA, hybrid functionals, and meta-GGA functionals). 

A new era in the computational chemistry of graphene/metal systems began with the de- 
velopment of approaches allowing to account for dispersion interactions, which are extremely 
important when investigating layered materials. A variety of methods were suggested to mit- 
igate the inherent problem of locality in DFT (for a review, see Ref. [H2])- Here we will name 
only approaches which have been employed so far for the investigation of graphene/metal 
interfaces (see e.g. Refs. [TTl 133]): (i) the dispersion correction schemes (referred to as 
DFT-D) by Grimme [31] and Tkatchenko and Schemer [35] ; (ii) utilization of truly nonlocal 
functionals, e.g. the so-called vdW-DF [36] . 

III. EXAMPLES OF THE GRAPHENE/METAL INTERFACES 
A. Graphene on Ni(lll), Ir(lll), and Rh(lll) 

As an example of the graphene/metal system let us consider the representative case of this 
interface, namely graphene/Rh(lll). In the unit cell of this system the (12 x 12) graphene 
layer is arranged on the (11 x 11) Rh(lll) stack [18j [191 ED EI]- This example we will use 
for the discussion of all features of the graphene-based interface with metal. Generally, this 
system is a good example of the lattice-mismatch interface between a graphene layer and 
the metallic surface. Its crystallographic structure is shown in Fig. [2} Considering the ar- 
rangement of the graphene layer on Rh(lll) one can identify several high-symmetry stacking 
positions for carbon atoms in the layer. (Different notations are used in the literature to 
mark these positions.) They are: 

• ATOP-position - carbon atoms surround the metal atom of the top layer and are 
placed in the hep and fee hollow positions of the Rh(lll) stack above (S — 1) and 
(S — 2) Rh-layers, respectively (hep — fee position) [circle in Fig. |2](a)]; 

• FCC-position - carbon atoms surround the fee hollow site of the Rh(lll) surface 
and are placed in the top and hep hollow positions of the Rh(lll) stack above (S) and 
(S — 1) Rh-layers, respectively (top — hep position) [rectangle in Fig. [2|a)]; 
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• if CP-position - carbon atoms surround the hep hollow site of the Rh(lll) surface 
and are placed in the top and fee hollow positions of the Rh(lll) stack above (S) and 
(S — 2) Rh-layers, respectively (top — fee position) [triangle in Fig. p[a)]; 

• BRIDGE-position - carbon atoms are bridged by the Rh atom in (S) layer [star in 
Fig-ga)]. 

All these arrangements of carbon atoms on transition-metal surfaces can be easily identified 
in the experiment, for lattice-matched as well as lattice-mismatched systems. The represen- 
tative examples are shown in Fig. |3j where the STM images of a graphene layer on Ni(lll), 
Rh(lll), and Ir(lll) are shown. The corresponding LEED images are presented as the 
insets. 

The graphene/Ni(lll) system is an intensively-studied representative example of the 
lattice-matched system. The STM and LEED experiments show the commensurate (1 x 1) 
structure of carbon atoms on Ni(lll) (Fig. |3j left) [21]. The lattice constant extracted 
from the experiment is 2.4 ± 0.1 A, which agrees well with the in-plane lattice constant of 
graphite (2.464 A). Earlier studies of this system considered the first intuitively prospective 
arrangement of carbon atoms, namely ATOP. This model was supported by available 
experimental studies of graphene/Ni(lll) by means of surface extended energy-loss fine- 
structure spectroscopy [38J. These experiments also found the relatively large distance of 
2.46 A between a graphene layer and Ni(lll) without any intermixing of graphene and Ni 
valence band states, that was in contradiction with available photoemission experiments. 
However, the later studies performed by means of LEED I(V)-cmve analysis found that 
HCP (top — fee) arrangement of carbon atoms above Ni(lll) is more realistic with the 
mean distance between carbon and top Ni layers of 2. 135 A [39]. This arrangement was 
later supported by many electronic structure calculations, ranging from LSDA to GGA and 
GGA-D and now this model is widely accepted for the arrangement of carbon atoms on 
the Ni(lll) surface [TTJ [T2J ESI 00]. The similar conclusions were made for graphene on 
Co(0001) p3] and Fe/Ni(lll) [TJJ HI]. 

The relatively small distance between graphene and Ni(lll) leads to the signifi- 
cant intermixing of the valence band states of both materials: the strong hybridiza- 
tion between graphene ir and Ni3d states was observed in the experiment [Fig. |4j upper 
panel] [HI HH ED SOI ME5\ and confirmed by theoretical results PH [TJ1 SOI 0HJ 07] . Here 
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FIG. 4: ARPES intensity maps along the V — M — K — T path in the Brillouin zone for 
graphene/Ni(lll) (upper panel) and graphene/Ir(lll) (lower panel). 

the electronic structure of free-standing graphene is drastically modified: for the occupied 
valence band states the graphene-derived a and tc states are shifted to higher binding energies 
by leV and 2.4 eV, respectively, due to the charge transfer from Ni states onto unoccupied 
states of graphene; the graphene Dirac cone is not preserved due to the strong intermix- 
ing of Ni and graphene valence band states and the large gap is opened around the K 
point of the graphene-derived Brillouin zone. The corresponding changes are also observed 
for the unoccupied valence band states. Following the above description, presently, the 
graphene/Ni(lll) is related to the system where graphene is "strongly" bonded to metallic 
substrate. The analogous description is also accepted for other lattice-matched systems, 
graphene/Co(0001) and graphene/Fe/Ni(lll). 
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The lattice-matched graphene/Cu(lll) system prepared via intercalation of thin Cu layer 
underneath graphene can be related to the "weakly" interacting system due to the completely 
filled d shell. Here the recent theoretical and experimental works show the weak interaction 
of graphene 7r and Cu valence band states j9j GU |43]. Also the Dirac cone is preserved in 
the system. The similar effects were observed for the system where graphene was grown by 
molecular-beam epitaxy of carbon atoms on single crystalline Cu(lll) |48j. 

In the case of lattice-mismatched graphene-based systems two cases can be identified: (i) 
"weakly" bonded graphene to metal, graphene/Ir(lll), graphene/Pt(lll) and (ii) "strongly" 
bonded graphene to metal, graphene/Ru(0001), graphene/Rh(lll). 

The graphene/Ir(lll) is a widely studied example of a "weakly" bonded graphene on 
the lattice-mismatched metallic substrate. Its representative STM image is shown in Fig. [3] 
(right image) in the inverted contrast imaging when ATOP positions are viewed as dark 
and FCC and HCP as a bright regions [49H51]. The conclusion regarding the "weak" 
type of interaction in the system under consideration was made on the basis of its almost 
intact electronic structure compared to free-standing graphene: (i) the binding energy of the 
graphene-derived 7r states at the T point is 8.15 eV, (ii) the 7r band has a linear dispersion 
around Ep, (iii) there is a small p-doping of graphene in this system [Fig. |1| lower panel] [521 
[53] . Additionally the moire superstructure of graphene on Ir(lll) results in modulation of 
the periodic atomic potential around carbon atoms followed by the appearing of the Dirac 
cone replicas and the opening of minigaps in the band structure of a graphene layer [52J. 
Experiments performed on this system by means of x-ray standing wave (XSW) technique 
yield a mean distance between graphene and Ir(lll) equal to 3.38 A with a corrugation 
between 0.6 and 1.0 A [17J. From the theoretical point of view the LDA approach gives the 
surprisingly good result for the mean distance between graphene and Ir(lll) of 3.42 A [541 
[55] . although one has to be aware that this can be due to the error cancellation and this 
method has a tendency to overestimate the binding. The GGA calculations, on the contrary, 
yield a very small binding energy of only several meV per carbon atom with a relatively 
large mean distance between graphene and Ir(lll) [4*9] 156]. The recent development of the 
methods which can take into account the non-local van der Waals interactions allowed to 
describe the graphene/Ir(lll) system correctly [T71 [57]. The obtained mean distance and 
the corrugation of the graphene layer are in the reasonable agreement with experimental 
data. In this DFT-D description the graphene layer is bonded to Ir(lll) via the van der 
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FIG. 5: Carbon atom-projected total density of states (a and tt) in the valence band for the differ- 
ent high-symmetry positions of the graphene/Rh(lll) system. The inset shows the corresponding 
density of states for the p z character only. 

Waals interaction with an antibonding average contribution from chemical interaction [T7] . 
These calculations indicate the charge accumulation at the graphene/Ir(lll) interface with 
the small charge transfer of ~ 0.01 electrons/C from graphene to Ir resulting in a slight 
p-doping of graphene with the Dirac cone shifted by 0.2 eV above Ep, which is consistent 
with experimental value of 0.1 eV fT7\ 152"] . The similar observations and conclusions are 
also valid for the "weakly" interacting graphene/Pt(lll) system [581 - I64"] . The later results 
are also consistent with the model calculations for the doping level taking into account the 
difference in the work function of graphene and the metal support [9l [26] . 

Graphene on Rh(lll) is chosen by US cLS cL typical example of the "strongly" bonded 
lattice-mismatched system. Its STM and LEED images are shown in the middle part of 
Fig. [3] [TSl [HJJ The moire structure is clearly visible in the experiment. The lattice con- 
stant of the moire structure extracted from the experimental data is ~ 29 A corresponding 
to the arrangement of (12 x 12) graphene layer on (11 x 11) Rh(lll). All high-symmetry ad- 
sorption sites of carbon atoms on Rh(lll) are marked on STM image by the corresponding 
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capital letters and these positions can be unambiguously assigned on the basis of compar- 
ison of calculated and experimental STM data [T9"] IBT)] . The measured corrugation of the 
graphene layer on Rh(lll) depends on the STM imaging conditions and was found to vary 
in the range 0.5 — 1.5 A depending on imaging conditions [HJ [19], [65] . The DFT-D optimized 
structure of graphene/ Rh(lll) is presented in Fig. [2|a) and the variation of the hight of 
the carbon atoms is shown in Fig. [2^b). Carbon atoms in the ATOP configuration define 
a high-lying region sitting at do = 3.15 A above Rh(lll), and those in other configurations 
form a lower region. The buckling in the graphene overlayer is 1.07 A. Carbon atoms in 
the BRIDGE configuration form the lowest topographic area 

(d = 2.08 A). The HCP 

and FCC regions are approximately 0.4 A and 0.8 A higher than the minima. Thus, the 
theoretically obtained value of corrugation is in very good agreement with those obtained 
in atomic force microscopy (AFM) measurements of this system of 1.1 A [19]. Authors of 
Ref. [19] have estimated the influence of dispersion forces on the obtained results: while 
qualitatively the observed picture remains the same, non inclusion of the van der Waals 
interactions (i.e. standard DFT-PBE treatment) yields larger corrugation 1.8 A) with 
a very similar low region (do = 2.10 A), but a high region at do = 3.90 A. This is due to 
the alternating "weak" and "strong" interactions of graphene with the Rh(lll) surface. In 
the case of the "strong" interaction between graphene and metal, standard GGA-treatment 
gives a reasonable result, whereas for the areas of the "weakly" bonded graphene dispersion 
forces, neglected by the standard procedure, are important. 

The coexistence of the "weakly" and "strongly" bonded graphene on Rh(lll) is reflected 
in the electronic structure. As can be seen from Fig. |5j where projected DOSs on the 
representative C atoms are shown, the low-lying region of graphene/ Rh(lll) displays rather 
a strong interaction involving hybridization of graphene states with those of the metal, 
whereas for the high-lying region a free-standing graphene-like picture can be observed 
(although small doping is visible). The difference in DOS curves around the Fermi level is 
not so pronounced as expected due to the fact that for the ATOP positions of a graphene 
layer it is close to the free-standing one with the low DOS at Ep (the finite DOS is due 
to the residual interaction with the Rh substrate). In case of the BRIDGE positions, the 
interaction is so strong that it might open a large gap in DOS for the graphene 7r-states. But 
the number of hybrid interface states can lead to the increasing of DOS in the vicinity of the 
Fermi level. This situation is analogous to the graphene/Ni or graphene/Co lattice-matched 
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(a) Band structures: graphene/Ni(111), spin-down 
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(c) Difference electron density plot: graphene/Rh(111) 
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FIG. 6: (a) Electronic band structure for different arrangements of a graphene layer on Ni(lll). 
(b) and (c) show the difference electron density, A / o(r) = p gr /Rh — /°Rh — Pgr, plotted in units of e/A 3 
calculated for different positions of graphene/Ni(lll) and graphene/Rh(lll). Red (blue) colors 
indicate regions where the electron density increases (decreases). 



systems, where such situation is observed [Fig. [6] (a)] [HI [12], [161 HO]- The difference in the 
DOS peak positions for the graphene a- and 7r-states reflects the strength of hybridization 
and the charge transfer in the system. Analogous to the above description, the strong 
interaction between graphene and Rh around the BRIDGE positions can be explained by 
the strong hybridization between graphene 7r and Rh d states in the valence band. That 
can lead to the partial charge transfer from Rh to graphene. This effect is reflected in 
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FIG. 7: Side and top view of the crystallographic structures of the (a) graphene/Fe(lll) and (b) 
graphene/Al/Ni(lll) intercalation-like systems. The corresponding units cells are marked in the 
figure. 

DOS as a shift of 7r-states of graphene from « 6.5 eV binding energy for ATOP position 
to higher binding energies for BRIDGE positions (several DOS peaks in the energy range 
E — Ep = —7.5... — 9.5 eV). The peaked resonances at the Fermi level can be viewed as the 
signature of an active dangling bond. This suggests that the HCP and FCC area display 
the strongest chemical activity in the moire overstructure and therefore one might expect 
preferable formation of TM clusters above these regions [37J. 

Now we will try to make a connection between electronic structure and geometry of 
lattice-matched and lattice-mismatched graphene/metal systems. For the graphene/Ni(lll) 
interface the electronic properties, strength of hybridization and level of doping of graphene 
depend on the arrangement of the system. For the top — fee [HCP) and top — hep [FCC) 
arrangements of carbon atoms on Ni(lll) the "strong" interaction and significant hybridiza- 
tion of graphene rr and Ni 3d states are observed [Fig. [6|a), left and middle plots]. Here, these 
effects lead to the opening of the energy gap between graphene ir and 7r* states with appear- 
ing of many, so called, interface states in this gap. For the fee — hep {ATOP) arrangement 
of carbon on Ni(lll) the Dirac cone is preserved and in this case the graphene layer is even 
slightly p-doped by 200 meV above Ep [Fig. |6](a), right plot]. Surprisingly, this value of 
doping is close to the one for the graphene layer on Ir(lll) [52J (see also discussion above 
in the text). We would like to emphasize, that for all possible arrangements of graphene on 
Ni(lll) the symmetry for two carbon sublattices in the graphene unit cell is broken, meaning 
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that the opening of the band gap in the graphene spectra is expected in all cases. However, 
as is evident from Fig. [6] the breaking of the symmetry is not the only prerequisite for the 
opening of the gap in the graphene layer. The strength of the additional modulation poten- 
tial (in this case, the strength of site-specific hybridization between valence band states) is 
important for the gap opening. This effect of hybridization, space and energy overlapping, 
of the graphene it and Ni 3d valence band states is clearly demonstrated in Fig. [6]^a,b) where 
for the top — fee (HCP) and top — hep (FCC) configurations the hybridization states ap- 
pear between a graphene layer and the Ni(lll) surface. For the fee — hep (ATOP) case 
the distance between graphen and Ni (111) is large preventing space overlapping of orbitals 
leading to the preservation of the Dirac cone. As an outlook, the two factors, the real space 
and the energy overlapping of the hybridized states defines the appearance and the value of 
the energy gap in graphene. 

For the graphene/Rh(lll) interface the situation is more complex. Here, all possible 
high-symmetry positions (see discussion above) are possible. In this case one can expect 
the different local hybridization strength as well as local electron doping of the graphene 
layer. This point is supported by the electronic structure calculations for this system [Figs. [5] 
and pj^c)]. The calculated local DOS for the graphene/Rh(lll) system is discussed earlier 
and it gives the idea about the local hybridization strength of the valence band states of 
graphene and Rh and the local doping. Considering the local charge distribution map 
[Fig. [6^c)], the similarity between systems, graphene /Rh(lll) and different arrangements in 
graphene/Ni(lll), is clearly visible. The ATOP, HCP, and FCC stackings of graphene 
on Rh(lll) and Ni(lll) have practically the same distribution of the electron density: (i) 
formation of interface states and (ii) charge transfer from metallic substrate to graphene for 
HCP and FCC, (iii) strong charge localization at the interface between graphene and metal, 
(iv) absence of any visible interaction for the ATOP positions. Therefore, in the case of 
graphene/Rh(lll) one may expect that the local electronic structure for these places will be 
similar to those of respective arrangements of graphene/Ni(lll) [Fig. |6](a)]. The most locally 
strongly interacting place for the graphene/Rh(lll) system is the BRIDGE position, which 
is, however, is not stable for graphene /Ni(lll) [33]. For this place the charge transfer and 
the interface charge localization is even stronger compared to the HCP and FCC places 
[Fig. [6|c)]. The nearly similar situation is realized for the graphene/Ru(0001) system, where 
the "strong" interaction between graphene and Ru is observed as well [22J. The opposite 
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case is graphene/Ir(lll) where interaction is "weak" and almost no hybridization between 
graphene and Ir valence band states as well as small charge transfer are observed [17] • 
However, we would like to emphasize that for all lattice-mismatched graphene /metal systems 
all above-listed high-symmetry positions for carbon atoms are realized. 

Summarizing all these facts one can expect that for the lattice-mismatched "strongly- 
interacting" systems in the electronic structure extracted from, for example, photoemission 
experiments all possible replicas of all crystallographic stackings have to be detected: two 
or more parabolas for graphene 7r or a states. However, available experimental data demon- 
strate the absence of such structures and single band always exist either for it or for o 
graphene-derived states [221 E21 E2] E7J [68]. These observations can be explained by the 
metallic nature of graphene in these systems. The high mobility of electrons and that 
during description of the electronic structure such lattice-mismatched systems have to be 
considered as a whole object without separation on "weakly" and "strongly" bonded regions 
lead to the fact that a graphene layer becomes fully doped via the charge transfer from the 
metallic substrate. 

The broken symmetry in the carbon lattice as well as locally distributed strong hy- 
bridization between graphene and substrate valence band states lead to the appearing of 
the relatively large gap around the K point of the graphene-derived Brillouin zone. This 
picture seems to be valid for "strongly" bonded graphene on Rh(lll) and Ru(0001), where 
in FCC and HCP regions the symmetry for two carbon sublattices is locally broken and 
strong hybridization of the graphene ir and d states of a substrate exists. For the "weakly" 
bonded graphene on Ir(lll) the local interaction between graphene and Ir substrate is small 
everywhere in the moire lattice and no strong hybridization is observed, meaning no gap 
around the K point. The slight doping of graphene in this system is defined only by the 
charge transfer from graphene on Ir(lll). The perturbation of the graphene-lattice poten- 
tial by the weak periodic potential of the moire lattice leads only to the appearing of weak 
replicas in the photoemission spectra with small gaps where they intersect with the main 
photoemission branches from graphene. 
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FIG. 8: Minority-spin electronic band structures calculated for the (a) graphene/Fe/Ni(lll) and 
(b) graphene/Al/Ni(lll) intercalation-like systems. The graphene 7r character in both plots is 
emphasized by thicker lines. 

B. Graphene-metal-based intercalation-like systems 

The properties of epitaxial graphene/substrate interfaces can be efficiently controlled by 
inserting other materials between graphene and the original support, a process referred to as 
intercalation. In case of the graphene /metal system, different examples were investigated, 
which can be found in the recent review papers [T3HT5] . Considering the graphene/Ni(lll) 
interface as a reference system, we would like to present here a few illustrative examples and 
the first one is dealing with the graphene/Fe/Ni(lll) intercalation-like system. The lattice 
structure, which was revealed by comparative analysis of the experimental (LEED) and 
theoretical (DFT/PBE) results [12], is shown in Fig. [7] (left): Fe atoms are placed in the fee 
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hollow sites following thereby the Ni(lll) structure and C atoms are staying in the top-fcc 
arrangement with respect to the Ni(lll) substrate. The interaction between graphene and 
metal is stronger here than in the graphene/Ni(lll) case due to the less filled d shell of Fe 
compared to Ni. Intercalation of iron underneath of the graphene layer changes drastically 
the magnetic response from graphene that is explained by the formation of the highly spin- 
polarized 3d z 2 quantum-well state in thin iron layer, keeping the electronic structure below 
the Fermi energy almost intact [compare Figs. [6] (a) and [8] (a) for theory (DFT/PBE), and 
Figs. [I] and [9j for experiment (ARPES)]. 

Thus, in the previous example, the unique conical electronic structure of graphene was 
completely destroyed. One can try to restore it via decoupling the graphene layer from metal, 
Ni(lll), by intercalation of Cu or Au, for example [211 HSl E9J [70] . Indeed, one may observe 
the desired picture, although the Dirac point is usually shifted up- or downwards with respect 
to the Fermi level, resulting in p- or n-doped graphene, respectively. Additionally, 7r states 
of graphene are hybridized with (i-states of the noble metal [2H ED] . The latter effect can 
be avoided, if necessary, by replacing Cu-Au with sp-metal, like Al. 

The crystallographic arrangement of atoms in the graphene/Al/Ni(lll) system is pre- 
sented in Fig. [7j Contrary to the graphene/Fe/Ni(lll) trilayer, both occupied and un- 
occupied states are modified by means of Al-intercalation [4TH [7T| [72] . As one can see 
from Figs. [4] and [8j in this case all electronic bands of graphene are shifted to lower bind- 
ing energies, compared to graphene/Ni(lll). Furthermore, the electronic structure of the 
graphene layer as well as the Dirac cone in the vicinity of Ep are fully restored (there 
is a small electron doping of graphene leading to a shift of the Dirac point below Ep 
by ca. 0.64 eV). The decoupling process in this case is also easily visible when comparing 
the near edge x-ray absorption fine structure (NEXAFS) spectra obtained for graphite, 
graphene/Ni(lll) and graphene/Al/Ni(lll) presented in Fig. [9j In the case of our reference 
system, graphene/Ni(lll), in the region of the Is — > n* transition, this spectrum has a 
double-peak structure compared to the one of graphite that is explained by the transitions 
of C Is core electron into two unoccupied states (interface states), which are the result of 
hybridization of the graphene and Ni valence band states [44, 45J. Additionally, one observes 
a reduction in the energy separation between the ir* and a* features compared to that in 
the spectra of graphite (a result of the lateral bond softening within the adsorbed graphene 
monolayer). Intercalation of thin Al layer underneath a graphene layer on Ni(lll) leads 
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to drastic changes in the C K NEXAFS spectrum: The shape of the spectrum, positions 
of main spectroscopic features as well as the energy separation between n* and a* features 
become similar to those in the spectrum of pure graphite. These facts immediately indicate 
that the graphene layer is decoupled from the substrate by intercalated Al. An interesting 
feature of the system under consideration is that owing to the symmetry breaking in the 
graphene lattice there is no energy gap for the ir states around the K point of the Brillouin 
zone of graphene and the both theory and experiment agree in this statement. At the same 
time, appearance of such a gap was recently demonstrated by ARPES for the graphene layer 
on Cu(lll), Ag(lll) and Au(lll) [221 I2H US]- Note, however, that these experiments are 
not supported by band-structure calculations where a very small [73] or no energy gap was 
observed for these surfaces [9]. 

Thus, by intercalation of different metals underneath graphene on Ni(lll) one can (i) 
shift the graphene unoccupied states away from Ep, when employing Co or Fe [T2l [71] or 
(ii) decouple the graphene states from the Ni(lll) substrate, when employing noble metals 
or Al [221 12H EH H21 EH] • In this way, in the latter case hybridization between graphene n 
and metal valence band states in the vicinity of Ep is completely avoided. 

IV. DISCUSSION 

Coming to the end of the present manuscript we will rise and discuss two main problems in 
the studying of graphene layers on metallic substrates. Going through the presently available 
literature which is mainly summarized in several review papers [131415] one can come to 
the conclusion that this field of graphene surface science research gives probably the only 
way for the preparation of high-quality graphene layers which can be either used in future 
devices in the contact with the supporting metal or can be transferred on the insulating or 
semiconducting support without losing the fundamental properties of graphene [3H7J. Here 
the controllable growth as well as high quality of obtained graphene layers open the door 
for the wide implementation of these technology in industry [751 EI| . Along with this rapid 
progress in the application of graphene, the nature of interaction of graphene with metals is 
far from the full understanding. 

As it was pointed earlier, initially the application of standard DFT methods (LSDA, 
GGA-PBE) gave in most cases the correct description of the geometrical arrangement of a 
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FIG. 9: (a) ARPES photoemission intensity maps for graphene/Fe/Ni(lll) and 

graphene/Al/Ni(lll) presented for the V — K direction of the graphene Brillouin zone, (b) C 
K edge NEXAFS spectra of graphene/Ni(lll), graphene/Al/Ni(lll), and a graphite single crys- 
tal. 

graphene layer on metallic substrates. However, the bonding strength was either overesti- 
mated (LSDA) or significantly underestimated up to non-bonding situations (GGA-PBE). 
When corrected for the missing dispersion interaction, the accuracy of DFT is currently be- 
coming acceptable and one can predict the bonding energies as well as the distances between 
graphene and metallic substrates (for lattice-matched and lattice-mismatched systems). 

The separation of the graphene/metal systems on two subclasses, lattice-matched and 
lattice-mismatched, does not give an answer on the different strength of interaction. In the 
subclass of lattice-mismatched systems both situations, "strongly" and "weakly" interact- 
ing graphene, are observed [graphene/Rh(lll) vs. graphene/Ir(lll)]. The first attempt 
to describe the difference in the bonding between graphene and metals was made in the 
framework of the so-called d band model |77| . which predicts a stronger binding with de- 
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creasing occupation of the d band giving a stronger binding when one goes in one group 
from the 5d to the 3d metals. According to this model the changes in the properties of the 
graphene/metal interface have to be more gradual as it is observed for the 3d and 4d metals. 
But the drastic changes in the bonding between graphene and Rh and Ir or graphene and 
Pd and Pt are observed. For example, the DFT studies by means of GGA-PBE predict that 
graphene is bonded to Ru, but repulsive interaction is predicted for adsorption of graphene 
on PtRu 2 , Pt 2 Ru, and Pt surfaces as well as on Ir J78J US]- This effect probably could be 
understood in terms of orbital overlap. In going from 4d to 5d elements, relativistic effects 
become significant, resulting in contracting the 6s orbitals and expanding the 5d orbitals. 
Probably, the overlap of the graphene 7r states with more diffuse 5d orbitals is worse than 
with more localized 4d orbitals. Later on, the effect of bonding of graphene on Pt(lll) 
and Ir(lll) was checked by the inclusion of van der Waals forces in calculations [UJ 180] . 
However, even these state-of-the-art calculations do not give a clear intuitive answer on the 
bonding mechanism between graphene and metallic surfaces. 

In this situation, the so-called incremental scheme |81| . which belongs to the group of 
local correlation methods [see e.g. Refs. [82| 183]]. may help considerably [MUST] . An asset 
of this approach is the possibility to analyze the individual correlation contributions to the 
binding energy, i. e. one can separate interactions between, e. g., s- or <i-orbitals of metal and 
those of graphene, that will give a better insight in the nature of bonding between graphene 
and substrate. 

Another interesting point concerning the metal/graphene interface is the existence of the 
energy gap at the K point in the graphene 7r-electron spectrum. The appearance of such 
gap is closely related to the violation of the symmetry in the graphene lattice. For the 
"strongly" interacting graphene/metal interfaces [graphene/Ni(lll), graphene/Co(0001), 
graphene/Ru(0001)] the existence of such gap is unequivocally described by the theory, 
both GGA-PBE and GGA-PBE-D, that also supported by the experimental observations. 
Here, as it was described earlier, the gap opening is due to the broken symmetry for two car- 
bon sublattices along with the strong hybridization between graphene n states and ci-states 
of the substrate. 

For the "weakly" interacting graphene/metal interfaces the situation is not so obvious. 
In the beginning we would like to recall the situation on the existence of such gap for the n- 
doped graphene on SiC(OOOl) [SHIES]- It was demonstrated that the appearance of such gap 
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in the spectrum of graphene, which was observed in Ref. [89] and assigned to the violation 
of the sublattice symmetry, can be assigned to the islanded graphene sample which could 
be obtained during underanneling of SiC substrate [BBJ • Moreover, the additional indication 
of the broken symmetry in the graphene layer is the observation of the intensity from the 
graphene 7r band in the second Brillouin zone along V — K direction. As was shown by 
Shirley et al. [90], if the potential for both carbon sublattices is identical, then due to the 
interference between emission from A and B sites of the graphene unit cell the intensity 
in the second Brillouin zone vanishes. In opposite case the final photoemission has to be 
detected. 

In the simple "weakly" interacting graphene/Ir(lll) and graphene/Pt(lll) the energy 
gap around K was not detected in experiment or predicted by theory j52j [61]. It was 
possible to open a gap only via strong hydrogenation of a graphene layer on Ir(lll) when 
local transformation from sp 2 to sp 3 hybridization of carbon occurs [91] or via heavy co- 
doping with Ir and Na [02]- In the later case this rehybridization of carbon atoms due to 
adsorption of Ir in HCP sites of moire structure can be locally enhanced by strong charge 
transfer in the presence of adsorbed Na. 

In case of the graphene/metal intercalation-like systems the earlier experiments demon- 
strate either existence of the gap, like for graphene/Na,K,Cs/Ni(lll) [93J, or its absence, 
like for graphene/Cu,Au/Ni(lll) 1 1 3 . 69]. However the correctness of these old results can 
not be taken as a reference due to the experimental limitations in former time (limited en- 
ergy and angular resolutions), because the fact of detection of such gap is very crucial to 
the method used. For example, later experiments performed with a modern display-type 
photoelectron detectors on graphene/K/Ni(lll) system demonstrate that intercalated K 
leads to heavy n-doping of graphene, decoupling of graphene from substrate, lack of any 
hybridization of graphene ir and substrate valence band states and the absence of any gap 
around the K point for the graphene 7r states [20J . 

The recent ARPES experiments performed on the graphene/Cu/Ni(lll), 
graphene/Ag/Ni(lll), and graphene/Au/Ni(lll) intercalation-like systems [21] indi- 
cate an energy gap around K of 180 meV, 320 meV, and meV, respectively. This is related 
by the authors of this work to the broken symmetry in the system, after noble metal 
intercalation, connected with the doping level of graphene in the obtained system. However, 
the available band-structure calculations on the level of LSDA do not reproduce these 
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observations [9j [26]. First of all this discrepancy between theory and experiment can be 
attributed to the non-correct description of the graphene/metal interface in the framework 
of LSDA which always overestimates the binding. But in this case the theoretically 
predicted splitting has to be even larger. Also, the recent work on the graphene/Al/Ni(lll) 
system do not show any gap for tt states around the K point in experimental as well as 
in theoretical data jlQ], although the doping level of graphene is quite significant with 
the Dirac point localized by ~ 0.64 eV below Ep. The existence of the energy gap for tt 
states was also demonstrated for graphene/Au/Ru(0001) system [23] that was supported 
by the observation of the non-vanishing intensity of the tt band in the second Brillouin 
zone (see discussion above). Nevertheless, this gap was detected after depositing on the 
initial graphene/Au/Ru(0001) system (initially is p-doped) K atoms, that was attributed 
in Ref. [21] to the effect of potassium doping of the system. Here, we would like to mention 
also several recent works where no or extremely small gap (if any) was observed for the 
"weakly" bonded graphene/metal systems [H21 IMH^T] . Results of these works point out 
that for the formation of the energy gap for graphene tt states the effect of broken symmetry 
in the graphene unit cell between two carbon sublattices has to be accompanied with the 
effect of strong modulation of the carbon periodic potential due to the hybridization of tt 
states with the valence band states of substrate. In the ideal case the space, energy, and 
wave vector conservations during hybridization have to be fulfilled. 



V. CONCLUSIONS 



The aim of the present perspectives article is to give a short overview of the main problems 
in the field of surface physics and chemstry of the graphene/metal interface. The main 
questions which are still open, concerning this system, are the origin of drastic variations of 
interactions in these systems (two subclasses) as well as an origin of modification of electronic 
structure of graphene in the vicinity of the K point. We have tried to organise these 
questions in a more ordered way and prepare roads for the solution of these problems. Several 
examples of the graphene/metal interfaces together with related intercalation-like systems 
were considered here from both experimental and theoretical points of view, that allows us 
to shed more light on the problem of the interaction of graphene with the substrates and 
modification of its electronic structure. At the same time, the recent progress in experimental 
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and theoretical investigations of the graphene/metal interfaces opens a wide door in the 
application of these systems in future electron and spin transport devices. 
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